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We show that the transfer of the plane wave spectrum of the pump beam to the fourth-order 
transverse spatial correlation function of the two-photon field generated by spontaneous parametric 
down-conversion leads to the conservation and entanglement of orbital angular momentum of light. 
By means of a simple experimental setup based on fourth-order (or two-photon) interferometry, 
we show that our theoretical model provides a good description for down-converted fields carrying 
orbital angular momentum. 
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I. INTRODUCTION 

It is well established that when the paraxial approxi- 
mation is vàlid, any electromagnètic beam with an az- 
imuthal phase dependence of the form e tl ^ carries an 
orbital angular momentum lh per photon 1]. Interest- 
ing enough by itself due to its fundamental character, 
this fact also raises possibilities for technical applica- 
tions. For example, in the rapidly developing field of 
quantum information, it has been pointed out recently 
that it is possible to increase the amount of information 
carried by a single photon byencoding qubits in the or- 
bital angular momentum 0,0. Laguerre-Gaussian (LG 
hereafter) beams are the most known and studied exam- 
ples of beams carrying orbital angular momentum. De- 
vices that discriminate the orbital angular momentum of 
Laguerre-Gaussian beams have been reported and exper- 
imentally tested for very low intensities, suggesting that 
they should work at the single photon level |4| . 

An important potential application of light beams car- 
rying orbital angular momentum is the generation of pho- 
ton pairs with discrete multidimensional entanglement 
0| • This can be obtained by means of spontaneous para- 
metric down-conversion (SPDC) pumped by a LG beam. 
Denoting by \m) a one-photon state carrying an orbital 
angular momentum mh and by / the azimuthal index of 
the the LG pump beam, the two-photon state generated 
by SPDC can be written as 



+oo 



guiar momentum, in most cases, the pump beam propa- 
gates in a birefringent nonlinear crystal as an extraordi- 
nary beam. The anisotropy of the medium causes a small 
astigmatism in the LG beam as it propagates, breaking 
its circular symmetry in the transverse plane. This sym- 
metry breaking is equivalent to an exchange of angular 
momentum between the medium and the pump beam, so 
that the conservation holds only on average. This effect 
depends on both the angular momentum of the pump 
beam and on the crystal length, being negligible for thin 
crystals and low vàlues of l. A detailed account of this 
problem will be published elsewhere. 

In an earlier paper we showed that the phase 
matching conditions in SPDC are responsible for a trans- 
fer of the amplitude and phase characteristics of the 
pump beam to the two-photon field. In fact, it is the 
plane wave spectrum, or the so-called angular spectrum 
of the pump beam that is transfered to the fourth-order 
spatial correlation properties of the down-converted field. 
In this work, we demonstrate theoretically and exper- 
imentally that the conservation of orbital angular mo- 
mentum as well as the multidimensional entanglement in 
the SPDC process in the thin crystal paraxial approxi- 
mation is a direct consequence of the transfer of the plane 
wave spectrum from the pump beam to the two-photon 
state. By means of a simple experimental setup based 
on fourth-order (or two-photon) interferometry, we show 
that our theoretical model provides a good description 
for down-converted fields carrying orbital angular mo- 
mentum. 
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This expression is based on the hypothesis that orbital 
angular momentum is conserved in SPDC. Some authors 
have studied this issue @, El 13 •> an d experimental re- 
sults 5] suggest that orbital angular momentum is in fact 
conserved. One has to consider, however, that although 
the process of down-conversion itself may conserve an- 



II. THEORY 

A. The state generated by SPDC 

In the monochromatic and paraxial approximations, 
the two-photon quantum sta te g enerated by non-collinear 
SPDC can be written as 0, [ÏÏf 
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where 

1^) = Yl C °s,cn JJ dq s dq_i $(q s , q») |q s , a s ) s \q u a i ) i . 

<Ta,<Ti D 

(3) 

The coefíicients C\ and C2 are such that | C 2 1 <C \C\\. 
C2 depends on the crystal length, the nonlinearity co- 
efíicient, the magnitude of the pump field, among other 
factors. The kets \c\j,<jj) represent one-photon states 
in plane wave modes labeled by the transverse compo- 
nent qj of the wave vector kj and by the polarization aj 
of the mode j = s or i. The polarization state of the 
down-converted photon pair is defined by the coefíicients 
C Gs ^ % . The function $(q s ,q^) is given by Q 



^(q s ,qi 



-q;) smc 



£|q 5 - guj 

4K 



, (4) 



where v(q) is the normalized angular spectrum of the 
pump beam, L is the length of the nonlinear crystal in 
the propagation (z) direction, and K is the magnitude 
of the pump field wave vector. The integrat ion domain 
D is, in principle, defined by the conditions q% < kg and 
Qi < kf. However, in most experimental conditions, the 
domain in which $(q s , q$) is appreciable is much smaller. 
If the crystal is thin enough, the sine function in Q can 
be approximated by 1. We assume that $(q s ,q^) does 
not depend on the polarizations of the down-converted 
photons. In some cases, this is not true, especially when 
one is dealing with type-II phase matching, in which case 
the two photons have orthogonal polarizations. However, 
this dependence can be made negligible by the use of 
compensators in the down-converted beams . 

The two-photon detection amplitude, which can be re- 
garded as a photonic wave function is 



*(r s ,r i ) = (vac|E t (+) (r i )E(+)(r s )|^), 



(5) 



where E^ + ^(r) is the field operator for the plane wave 
mode j. In the paraxial approximation, E^ (r) is 

^+)(r) = e ikz [ dqa j (q,cr)e a e i(q ·' , -é z ). (6) 
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The operator a j (q, <j) annihilates a photon in mode j 
with transverse wave vector q and polarization a. 

In the analysis that follows, we do not need to consider 
polarization. So, ^(r s ,r^) will be treated as a sealar 
function. In addition, we will work in the far field and 
make the following simplificat ions: z s = z\ = Z, k s = 
ki = \K. It is known that if the paraxial approximation 
is vàlid, the two-photon wave function is 

9{ PaiPi , ZB , Zi ) = s[^±Bi,z s j r(?^,z 

' (7) 

where £(p, z) is the normalized elèctric field amplitude 
of the pump beam and T{p,z) = sinc(|^p 2 ). 



In order to clean-up the notation, we will omit the de- 
pendence on the z coordinate hereafter. We see that 
the two-photon wave function \I/ carries the same func- 
tional form as the pump beam amplitude, calculated in 
the coordinate p = \p s + \pi- The pump beam field 
£(p) is characterized by its wavelength À G and its waist 
w Q . To be more precise, we will write £ as £ (p; \ ,w ). 
Since we are working with down-converted fields satisfy- 
ing X s = Xi = 2À G , it is convenient to write ^ in terms of 
a beam with the same angular spectrum of the pump 
field, as required by Eq. {7J), but with a wavelength 
A c = 2À , and a waist w c = \[2w . From the general 
form of gaussian beams, apart from a phase factor and 
normalization constants, it is evident that 



£{p\ \ 



£(V2p; 2A , V2w ) = U(V2p). (8) 



So, \I/ can be put in the more convenient form 



V2 
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Let us now suppose that the down-converter is pumped 
by a LG beam whose orbital angular momentum is l h per 
photon, described by the amplitude £ l p {p\ À , w Q ). Here, 
p is the radial index. In order to study the conservation 
of angular momentum in SPDC, we will expand the two- 
photon wave function \I/(p s , pj in terms of the LG basis 
functions U l p s s (pj U l p \ (pj. That is, 

*(p a , ft ) = EE^:W^^) ( 10 ) 

l s ,p s li, Pi 

From the orthogonality of the LG basis, C p s s 1 ^. is given by 



Cl't = lldp.dp, Ui ( Bl+B. )jríPs-P 
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xu; l ;{ Ps )u; l ;{ Pi ). (li) 

Let us make the following coordinate transformation 
in Eq. dl: R = p s + p- and S = \{p s - Pi ). So, 

C l p ± =JJdRdSU l p (J^j F(V2S) 

x u; 1 ; (| + s) u* - s) • (12) 

When L is small enough (the thin crystal approxima- 
tion), T can be approximated by 1 in Eq. (|T2|) . provided 
the order of the LG modes (N = 2p + |/|) is not too 
large. In this case, the integral in S is proportional to 
U p l s s (R) *Z/Ç^(R), that is, the convolution of U p [ s and 
Up* 1 .*. Numerical calculat ions show that in the worst case, 
that is, p s = Pi, \l s \ — \h\ an d R = 0, for a 1 mm-thick 
crystal, pumped by a làser with A = 351 nm and a waist 
of wo = 1 mm, the mean square error is less than 1% for 
N = 100. Since we are neglecting the effects due to the 
anisotropy of the crystal, as discussed before, there is no 
point in seeking exact solutions for large vàlues of N. 
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Under the thin crystal approximation, Eq. (|12|) is more 
convenient ly written in terms of Fourier transforms, as 

C l p ± cx JJJdRdqd^ v£(V2q') 

xV p t(q)V;f(q) e m <«'-^ 
oc ^q V l p (V2 q) (q) V p f (q), (13) 

where V^(q) is the Fourier transform of U^(H). Writing 
Eq. (|13|) in cylindrical coordinates q — > (q,(j>), the LG 
profiles are V^(q) = v^(q) e Lh> ^ . Then, we have 

Cp'X « JJqdqd<Pv l p (V2q)v; l ;(q)v^(q) e -*C+'«-'>* 

(14) 

that is, 

oc ^(>/2 9 )<-(<ï) (15) 

Thus, orbital angular momentum is conserved in the 
SPDC process. In principle, this conservation could be 
satisfied by fields exhibiting either a classical or quan- 
tum correlation (entanglement) of orbital angular mo- 
mentum. We will now show that the conservation leads to 
entanglement of orbital angular momentum of the down- 
converted fields. 

From (|ÜJ it is evident that, when T = 1, the bipho- 
ton wave function reproduces the pump beam transverse 
profile. Let us assume that (|ÜJ) (with T = 1) accurately 
describes the two-photon state from SPDC and that the 
pump beam is a LG mode with 1^0. Then, the biphoton 
wave function is 

*(p a ,p.)=Z4(£2±£) (16) 

from which, it is evident that ^(p s -f A,p^ — A) = 
^f(p s ,p i ). Due to the phase structure of U^, for l ^ 
there exist transverse spatial positions p s0 and p i0 such 
that Víp(p s o + Píq) = 0. Then, clearly 

V(p s0 + A, p i0 - A) = ^(p s0 , p i0 ) = (17) 

and the coincidence detection probability V(p s1 p i ) = 
\^{PsiPi)? satisfies 

V(p s0 + A, p l0 - A) = V( Ps0 , p l0 ) = (18) 

Now suppose that the down-converted fields exhibit a 
classical correlation that conserves orbital angular mo- 
mentum. The detection probability V cc for such a corre- 
lation can be written as 

oo 

r cc ( Ps , Pi )= Yl W-iM\ 2 \GiM\ 2 (19) 

h = — oo 

where Fi s (p s ) and Gi i (p i ) represent down-converted sig- 
nal and idler fields with orbital angular momentum l s h 



and lih per photon, respectively. Here the coefficients P\ i 

satisfy J2u=-oo P li = 1 and P h ^ °- Now > if 03 accu - 
rately describes the two-photon state, then it must also 
satisfy the equivalent of (|T%|) : 

V cc (p s0 + A, p i0 - A) = V cc (p s0 , p i0 ) = 

which gives 

oo 

P il |F ; _ íi (p s0 + A)| 2 |G íi (p i0 -A)| 2 = (20) 

li — — OO 

Since Pi. > 0, a non-trivial solution to (|2Uj) exists (for 
the cases where at least one P\ i ^ 0) only if \Fi-i.(p 30 + 
A)| 2 = or \G h (p i0 - A)| 2 = for all A, which implies 
that Fi-i i = or G^ = 0. Thus, a classical correlation of 
orbital angular momentum states cannot reproduce the 
two photon wave function 

With the reasoning above, we have shown that, as- 
suming Eq. (|ÜJ) accurately describes the biphoton wave 
function from SPDC, the conservation of orbital angular 
momentum in SPDC is not satisfied by a classical cor- 
relation of the down-converted fields. This implies that 
the fields are entangled in orbital angular momentum. 

B. The Hong-Ou-Mandel interferometer 

Having demonstrated that the two-photon wave func- 
tion J7J) leads to conservation and entanglement of orbital 
angular momentum, the next step is to prové that it de- 
scribes accurately the state generated by SPDC within 
the assumed approximations. Although direct coinci- 
dence detection provides information about the modu- 
lus of ^(p s ,pj, its phase structure can only be revealed 
by some sort of interference measurement. We do this 
with the help of the Hong-Ou-Mandel (HOM) interfer- 
ometer [T^ , represented in Fig. Q]and described below. 
Coincidence measurements are taken from the two out- 
put ports of the beam splitter. When the interferometer 
is balanced, that is, when paths s and i are equal, we 
have fourth-order interference. When the path length 
difference is much greater than the coherence length of 
the down-converted fields, the interferometer plays es- 
sentially no role other than decreasing the coincidence 
counts by a factor of 1/2, and we can perform simple 
coincidence measurements. 

In the HOM interferometer, the state J3J) is incident 
on a symmetric beam splitter as shown in Fig. The 
annihilation operators in modes 1 and 2 after the beam 
splitter can be expressed in terms of the operators in 
modes s and i: 

ai(q) =ta s (q x ,q y )+irai(q x ,-q y ) (21) 
a 2 (q) = tai(q x , q y ) + ira s (q x , -q y ), (22) 

where t and r are the transmission and reflection coeffi- 
cients of the beam splitter. The negative sign that ap- 
pears in q y components is due to the reflection from the 
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FIG. 1: The Hong-Ou-Mandel interferometer. Reflection at 
the beam splitter causes a sign change in the ?/-coor dinat e. 



beam splitter, as sl·iown in Fig. If ri and r 2 are the 
positions of detectors D\ and D 2 , each located at one 
output of the beam splitter, the coincidence detection 
amplitude is given by 



^ c = ^ u (ri,r 2 ) + ^ rr (ri,r 2 ), 



(23) 



where the indices tt and rr refer to the cases when both 
photons are transmitted or reflected by the beam splitter, 
respectively. Combining Eq. ® with T = 1 and Eqs. 
(|2T1-I25j). it is straightforward to show that, for t = r = 
1/a/2, apart from a common factor, 



# c (Pi,p 2 ) °<\ 



U 
U 



xi + x 2 yi + y2 
-2/1 - V2 



V2 

xi + x 2 



(24) 



V2 ' V2 

Since the pump beam is a LG beam, U has the form 

U{p)=u l p {p)e a +. (25) 

According to Eq. (J2U), the corresponding coincidence 
detection amplitude is 

* c (Pi,p 2 ) = * C (R,0) oc u l p (R)smW, (26) 

where R — ^|Pi + P 2 | and # is defined by the relations 

pi sin (/>i + p 2 sin 2 



sine/ = 



cos # 



i? 

pl COS 01 + P2 COS (j)2 

R 



(27) 
(28) 



The coincidence detection probability, which is propor- 
tional to |^ c (i?,#)| 2 , is 



Pi 2 (p l5 p 2 ) oc \u l AR)\ 2 sin 2 16. 



III. EXPERIMENT 



(29) 



The experimental setup, shown in Fig. consists of 
two bàsic parts. The first part is the generation of a 
Laguerre-Gaussian (LG) mode using a mode converter, 
which transforms a Hermite-Gaussian (HG) mode into a 
LG mode. A detailed account of mode conversion can be 
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FIG. 2: Experimental setup. A wire is inserted into the làser 
cavity in order to generate a HG mode. A mode converter 
consisting of two identical spherical (L s ) and cylindrical (L c ) 
lenses converts the HG mode to a LG mode of the same or- 
der. The LG mode is then used to pump the BBO crystal, 
generating entangled photons incident on a HOM interferom- 
eter. The beam splitter BS is mounted on a motorized stage. 
Coincidence counts are recorded at detectors D\ and D2. 



found in refs. E^H- To create the HG-mode, we insert 
a 25 /im diameter wire into the cavity of an Argon làser, 
operating at ~ 30 mW with wavelength 351.1 nm. The 
wire breaks the circular symmetry of the làser cavity. It 
is aligned in the horizontal or vertical and mounted on 
an £2/-translation stage. By adjusting the position and 
orientation of the wire, we can generate the modes HG01 , 
HG10, HG02, and HG 2 o- The beam then passes through 
a mode converter consisting of two spherical lenses (L s ) 
with focal length f s = 500 mm and two cylindrical lenses 
(L c ) with focal length fc = 50.2 mm. The first spherical 
Iens is used for mode-matching and is located « 1.88 m 
from the beam center of curvature. The second spheri- 
cal Iens is placed confocal with the first, and is used to 
"collimate" the beam. The cylindrical lenses are placed 
(d = f c / V% ~ 35 mm) on either side of the focal point of 
Iens L s and aligned at 45°. The cylindrical lenses trans- 
form the HG mode into a LG mode of the same order by 
introducing a relative tt/2 phase between successive HG 
components (in the ±45° basis, due to the orientation of 
the cylindrical lenses) of the input beam 0? Q • The 
quality of the output mode was checked by visual exam- 
ination of the intensity profile [l5( as well as by inter- 
ference techniques: using additional beam splitter s and 
mirrors (not shown), the interference of the LG pump 
beam with a plane wave resulted in the usual spiral in- 
terference pattern [Ï4| . 

The second part of the setup is a typical HOM in- 
terferometer 12]. The Argon làser is used to pump 
a 7 mm long BBO (/3-BaB 2 04) crystal cut for type II 
phase matching, generating non-collinear entangled pho- 
tons by SPDC. The down-converted photons are reflected 
through a system of mirrors and incident on a beam 
splitter with measured transmittance T « 0.67 and re- 
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íiectance R « 0.33. Since the down-converted photons 
are orthogonally polarized, a half-wave plate (HWP) is 
used to rotate the polarization of one of the photons 
(V — > H). A computer-controlled stepper motor is used 
to adjust the position of the beam splitter. The detec- 
tors are EG & G SPCM 200 photodetectors, mounted 
on precision translation stages. D2 remained fixed while 
a computer-controlled stepper motors were used to scan 
detector Di in the transverse plane. Coincidence and sin- 
gle counts were registered using a personal computer. In- 
terference filters (1 nm FWHM centered at 702 nm) and 
2 mm circular apertures where used to align the HOM in- 
terferometer. The transverse intensity proflles were mea- 
sured with the interference filters removed and circular 
apertures with diameter 0.5 mm and 1 mm on D\ and 
D2, respectively. 



IV. RESULTS AND DISCUSSION 

The results are shown in Figs. |3|to[Sl The left sides 
of the figures show the expected coincidence patterns, 
obtained from the squared modulus of Eq. JÜJ) in the 
non-interfering regime (interferometer unbalanced), and 
from in the fourth-order interference regime (inter- 
ferometer balanced). The right sides of the figures show 
the measured coincidences. 

In Fig. [3J the nonlinear crystal was pumped by a LGj 
(1 = 1) beam. Its intensity profile is shown in Fig. 
in agreement with Eq. @. In the interference regime, 
shown in Fig. Ob, the two interference peaks predicted 
by Eq. ((29)) are easily seen. In Fig. 21 the nonlinear 
crystal was pumped by a LGq (l = 2) beam. Now, the 
interference pattern shows four peaks, in agreement with 
Eq. (E||). 

In order to test the the translational invariance of 
*ff(p 1 ,p 2 ), which leads to the conclusion that the two- 
photon state is entangled in orbital angular momentum, 
we repeated the measurement of Fig. [3d, with detector 
D2 displaced by Ax = Ay = 1 mm. The interference 
pattern obtained is shown in Fig. [SJ The coincidence 
pattern measured by scanning D\ is now dislocated by 
Ax = Ay = — 1 mm, still in agreement with Eq. (|2ÜJ). 



V. CONCLUSION 

We have shown experimentally that our theoretical de- 
scription of the two-photon wave function is accurate. 
Information about its modulus and phase structure were 
obtained by direct coincidence detection and coincidence 
detection of fourth-order HOM interference effects, re- 
spectively. The transfer of the plane wave spectrum of 
the pump beam to the fourth-order transverse spatial 
correlation function of the two-photon field generated by 
SPDC leads to the conservation and entanglement of the 
orbital angular momentum of the down-converted fields. 
We should stress that this effect is restricted to the con- 
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FIG. 3: Coincidence proflles predicted (left) and measured 
(right) when the crystal is pumped by a LGj beam. a) No- 
interference regime (Hong-Ou-Mandel interferometer unbal- 
anced). b) Fourth-order interference regime (interferometer 
balanced) . 



text of two approximations. The first is the par axial ap- 
proximation, in which our model for the transfer of plane 
wave spectrum in SPDC is based. However, the paraxial 
approximation is also the context in which the angular 
momentum carried by electromagnètic beams can be sep- 
arated into an intrinsic part, associated to polarization, 
and an orbital part, associated to the transverse phase 
structure of the beam. The second approximation is the 
so-called thin crystal approximation. It is possible to 
show that this approximation would not be necessary if 
the non-linear medium were isotropic. The birefringence 
of the crystals used for SPDC causes non-conservations 
of the orbital angular momentum that are proportional 
to the crystal length. Rigorously speaking, orbital angu- 
lar momentum would never be conserved in SPDC due 
to this effect. In thin crystals (few millimeters in length) 
however, it can be neglected. We believe that the argu- 
ments and the experiment reported here provide addi- 
tional evidence of conservation and entanglement of the 
orbital angular momentum of light in SPDC, as well as 
the limits within which they should be understood. 
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FIG. 4: Coincidence profiles predicted (left) and measured 
(right) when the crystal is pumped by a LGq beam. a) No- 
interference regime (Hong-Ou-Mandel interferometer unbal- 
anced). b) Fourth-order interference regime (interferometer 
balanced) . 
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FIG. 5: Coincidence profile predicted (left) and measured 
(right) when the crystal is pumped by a LGj beam, in the 
fourth-order interference regime (Hong-Ou-Mandel interfer- 
ometer balanced) . Detector D2 was displaced by Ax = Ay = 
1 mm. 
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